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Abstract 

We use the chiral Lagrangian to investigate the global properties of the A-flavor QCD 
vacuum as a function of the 6 parameter. In the case of exact quark degeneracy we 
find evidence for first order phase transitions at 6 = n ■ (odd integer). The first order 
transitions are smoothed by quark mass splittings, although interesting effects remain 
at realistic quark masses. We emphasize the role of the ry' condensate in our analysis. 
Finally, we discuss the implications of our results on the internal hadronic structure of 
axion domain walls and axion cosmology. 
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1 Introduction 

The QCD Lagrangian may be augmented with the gauge invariant term 9 GG/327r^ which, 
despite being a total divergence, can have physical effects due to gauge field configurations 
with non-trivial topology. In nature this term appears to be absent (or at least small - 
9 < 10~^° [0), but it is of theoretical interest to determine the effect such a term has on 
the theory. It is possible that understanding the effects of a non-zero 9 angle may help 
determine why the term is absent in the strong interactions. In this paper we shall study the 
effects of adding a non-zero 9 angle on the QCD vacuum assuming it does not change the 
pattern of chiral symmetry breaking {SU{N)l ® SU{N)ji SU{N)y) in the quark sector. 
While it is possible that adding a non-zero theta term to QCD may have more dramatic 
effects, such as altering the nature of confinement @, ^, we consider ours to be the more 
conservative assumption. In fact, it can be shown |Q that some exactly solved softly broken 
supersymmetric gauge theories which confine by dyon condensation continue to do so for all 
9 []. If the behaviour of QCD is smooth with changing 9 then the low energy theory may be 
described in terms of the standard chiral Lagrangian formalism. Witten |p has previously 
studied this possibility and shown that there may be first order phase transitions at = tt. 
He was particularly interested in the possibility of spontaneous CP violation aX 9 = -k (note 
that although ^ 7^ in general violates CP, 9 = tt actually conserves CP). Creutz further 
explored the phase structure of the mass degenerate chiral Lagrangian, demonstrating the 
existence of first order phase transitions at ^ = nvr (n e odd integers). Here we develop these 
observations with further study of these phase transitions, concentrating in particular on the 
role of the r]' meson, and study their behavior in the presence of mass splitting in the quark 
sector. We find that even for realistic values of the quark masses interesting effects persist 
as 9 is varied from to 27r. Finally, we discuss the implications of these phase transitions on 
axion models, in which ^ is a dynamical variable. In particular, we investigate the internal 
hadronic structure of axion domain walls (previously studied by Huang and Sikivie 0), and 
their coupling to thermal degrees of freedom in the early universe. We conclude that the 
domain walls can efficiently dissipate their kinetic energy and therefore contract away very 
rapidly after the QCD phase transition. 

2 N Degenerate Quark Flavors 

As a first example we consider QCD with quark flavors and a non-zero 9 angle 

C = -^G^'^G,, + ^^it^D^ - M),,vi>^- + -l^e^,„^G"^^G-^ (1) 

^However, in these models the low energy effective 9 angle depends in a complicated fashion on the bare 
9 angle. 
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where i runs over the flavors, and we shall take M = ml. The bare Lagrangian has 
an SU{N)l ® SU{N)fi U{l)v global flavor symmetry weakly broken by the small current 
quark masses. The axial U(l) symmetry of the classical theory is anomalous in the quantum 
theory. Using the chiral Ward Identity associated with the U{1)a symmetry we may rotate 
the 9 term on to the quark mass matrix which becomes Mg = Me*^/^. 

We shall assume that for all 6 the theory conflnes and breaks the chiral symmetry of the 
quarks to the vector subgroup. At low energies the theory may then be described by the 
standard SU{N)l ® SU{N)ji chiral Lagrangian 

C = ^tr{df,U^d^'U) + ERetr{MeU^) , (2) 

where U{x) = exp{2m'^T°' / F) is an SU{N) matrix. Taking a derivative with respect to M^- 
in (|1|) and (0) we note that 

{^,^,) = , (3) 

so U can be directly related to condensates of quark bilinears. 

The behavior of the minimum of the potential with changing 6 is easily determined when 
N = 2 since the trace of a general SU (2) matrix is real and the potential reduces to 

eo = -m^ cos{e/2) RetrU = - cos{e/2) u , (4) 

where we have deflned u = RetrU. The variable u varies between —1 and 1, and for 
cos{6/2) 7^ the potential energy is minimized at one of those values. These vacua corre- 
spond to U matrices 

U-±(l']. (5) 



When 6 = (odd integer) ■ vr the potential is exactly zero (at lowest order in the potential) and 
all ?7's are exactly degenerate. It is easy to see that at these values of 6 there is a flrst order 
phase transition with the vacuum jumping discontinuously between the two conflgurations 
(^). We show the form of the potential V graphically in Figure 1 for varying 6 and for the 
set of conflgurations U = diag(e*", e~*") (non-zero a corresponds to a 7r° condensate). The 
vertical axis in the flgure is measured in units of Em. The flrst order phase transition at 
6' = TT is clear. 

Creutz [0 has investigated the global behavior for general N. In this case there are N 
distinct minima of the form 

^_g-2n«/ivj n = i^2,...N , (6) 
occurring for the special values of 9 

9n = 2n7r n = l,2,...N . (7) 
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These minima are connected by N first order phase transitions at ^ = (odd integer) it. We 
note that the minima correspond to the vector flavor symmetry preserving diagonal SU(N) 
matrices. The U{1)a symmetry provides a simple understanding of the origin of the vacua 
in (^), which are characterized by bilinear quark condensates. The U{1)a symmetry of (|l|) 
is broken by the anomaly to a Z2N discrete symmetry acting on the quark flelds. Therefore, 
discrete Zjv transformations on the quark condensates relate equivalent vacua of the theory. 




Figurel: Potential, scaled by Sm, for U = diag(e*°, e and varying 6 in the degenerate 

two flavor case. 



There is an important subtlety to discuss at this point. To include the effects of the 6 
angle in the chiral Lagrangian we have made a U{1)a rotation on the quark fields and moved 
the 6 dependence into the quark mass matrix. To put the low energy Lagrangians of bare 
theories with different 6 angles into the form of (0), different axial rotations have been made 
in each case. Before comparing the quark condensates in different models we must undo (or 
compensate for) this axial rotation. We denote the matrix of condensates in the original 
basis of (|I|) as Uq = e~*^/^f/. In the original basis the overall phase of the quark condensate 
changes smoothly with 6 and a discontinuous jump in the vr^ condensate occurs at ^ = vr. 
This is shown in Figure 2 for the case of = 2, where the Uq is parametrized as 

Uo = e-'^l^ { °. \ . (8) 

After including the compensatory axial rotation, the vacua at ^ = 0, 27r are now identical, 
and given by 
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Similarly, in the flavor case the compensatory rotations also cancel the phase factors in 
(^). Thus in the original quark basis the vacua found by Creutz are actually the same 
vacuum and correspond to Uq = X. 
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Figure 2: Phases of the vacuum Uo matrix as a function of 6 in the degenerate two flavor 

case. 

In the two flavor case, for = tt ± e, the (compensated) minima are respectively 



Since they have an associated phase they violate CP. This model realizes the possibility 
discussed by Witten that at = tt there may be spontaneous CP violation by the occurence 
of two CP violating vacua. (CP is not explicitly violated when 6 = tt since physics is invariant 
to 2tt shifts in 6, and CP acts by changing the sign of any phase in the theory, 6 = ±7r are 
equivalent.) 



To better understand the physical implications of the vacuum dependence on 6 it is helpful 
to consider the problem in an SU{Nc) gauge theory at large Nc. In this context the phase 
of U can be identifled with a condensate in the t]' fleld. It is well known that the anomaly 
in the U{1)a symmetry is the result of the quark triangle graph which is suppressed at 
leading order in l/Nc- Thus at large Nc the full U{N)l ® U{N)ji global symmetry on the 




(10) 



3 Large Nc 



4 



fermions is restored. The fermion condensates break the symmetry to U{N)v and there are 
A^^ goldstone bosons with the new state being identified with the t]'. At 0{1/Nc) in the Nc 
expansion the rj' acquires a small mass. If for a given large value of Nc we take the quark 
masses to zero then this subleading term can dominate the r]' mass. In |^ it was argued that 
this is the limit most appropriate to realistic QCD. We can hope therefore to obtain some 
insight into the transition described above by looking to the large limit. 
The effective theory at large Nc is p 

£ = ^tr^d^U^d^'U) + T.Retr{MU^) - ^{(j)-ef. (11) 

Here U contains the A^^ goldstone degrees of freedom, and is no longer restricted to unit 
determinant. The additional phase degree of freedom correponds to the r]' meson, with 



detU{x) = e-''^(") 

and 

\/2N 

m = -^v' ■ (12) 

To obtain (|lT]) we have performed an axial rotation on the bare Lagrangian to move the 6 
dependence into the mass matrix, then undone the rotation by absorbing 6 with a shift in the 
Tj' field. This redefinition of the phase of U corresponds exactly to the previously discussed 
subtlety of ensuring that we compare theories with different 6 in the same quark basis. It is 
important to note that the final term in ([lID is only the first term in an expansion in (p — 6 
of a function that is 27i periodic. 

Restricting our attention again to the case N = 2 the potential of the large A"c theory is 

V = -ERetr{MU^) +^{(j)-ey . (13) 

We are interested in the case rfEm — *■ oo, where the rj' is very heavy, and we expect to 
recover the analysis of the SU{2) chiral Lagrangian above. The second term in the potential 
dominates and is minimized by = 6'. As ^ changes between and 27r the phase of the mass 
matrix, which is the t]' condensate, behaves as in Figure 2. An ambiguity in the description 
arises here because a vr*^ condensate corresponding to a = vr has the same effect as an rj' 
condensate of = 27r on the matrix U in (P). This ambiguity is resolved in the presence 
of mass splittings, as the discontinuities are smoothed, and one can resolve the different 
behaviors of the vr*^ and rj' condensates. We discuss the effect of mass splittings below. 
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4 Mass Splitting 



Returning to the SU{N) chiral Lagrangian we may now consider the effects of sphtting the 
quark masses. The potential is then 



N-l N~l 

V = -^J2^i - 0/N) - Smjv cos(^/A^ + ^ «i 



(14) 



As before we have made an equal axial rotation on each quark flavor to shift the 6 angle 
onto the quark mass matrix. The minima are of the form 



{U) =diag[e"^\e 



(15) 



If we take nii^f/mi oo the last term in (|T^) dominates and the heavy quark condensate 



acquires a phase J2 



N-l 



9/N . In other words the phase of the A^th diagonal element of 



UMq goes to zero corresponding to decoupling of the heavy quark and the problem reduces 
to that in which the 9 angle has been rotated onto just the — 1 light quarks. This is why 
the strong CP problem cannot be solved by simply rotating the 6 angle solely onto a heavy 
quark such as the top. In nature the charm, bottom and top quarks are sufficiently heavy 
that they essentially decouple from the 6 dependence of the vacuum. 

For intermediate mass splittings we must minimize (p!4D. For N = 2 this is most easily 
done numerically. Minimizing ([I^) and rotating the resulting U back to the quark basis of 
(P) we have 

' 



Un 



'ie/2 







(16) 



We show the phase a that minimizes the potential at varying 6 for increasing mass splitting 
in Figure 3. The first order phase transition is smoothed out by the mass splitting and as 
fTT-d/fnu ^ oo we recover the one flavor case for which there is a single vacuum. 
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Figure 3: The phase a vs 6* in the presence of mass splitting 
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With mass splittings there is only a single CP conserving minimum at ^ = vr: 

This minimum violates the vector symmetry of the model, however, this is just a reflection 
of the explicit breaking of that symmetry in the non- degeneracy of the quark masses. Fur- 
thermore, in this case one can easily see by examining the form of U that both the t]' and 
the tt" condensates are changing with 6. 

The three flavor case is more complicated. The minimization equations reduce to (with 

m„sin(a — 9/3) = rrid sin(/5 — 9/3) = — nis sin(a + /? + 9/3) (18) 

The behavior is most easily studied by displaying the solutions at 6' = vr, where phase 
transitions are likely to occur. There are several types of solutions to (plsl). Trivial solutions 
occur when the sine functions are all simultaneously zero. For realistic quark masses (m„ ~ 
5 MeV, rrid ~ 10 MeV, ~ 150 MeV) one of these solutions is the global (and only) 
minimum corresponding to the heavy strange quark having decoupled and the minimum 
closely resembles the N = 2 case with ^ irid (See Figure 3). Non-trivial solutions were 
studied previously by Witten |^, and occur when the quark masses satisfy 

TTiumd > ms\md - mj\ , (19) 

which includes the degenerate mass limit discussed above. They satisfy 



1 - cos(a + /5 + ./3) = «-^^-^^--^uf (20) 



These latter solutions exhibit spontaneous CP violation. In this case the trivial solutions 
discussed above become points of inflection. Witten's solutions are probably excluded for 
realistic values of the light quark masses. 

As an explicit example we have examined the behavior with varying 9 of the Large- 
A^c effective Lagrangian. We used the physical 77,77' masses to determine r/S = 200 MeV 
and rris = 150 MeV. With such large values of r and the r and dependence of 
the solutions effectively decouples, and we are left with the behavior of the SU{2) chiral 
Lagrangian. For example, taking m„ = = 7 MeV, we recover the flrst order transition 
of the degenerate two flavor case. Minimization in the SU (3) degrees of freedom m. U = 
g-j0/3^^^^j^gj(o+/3)^ g-j(2/3)^ leads to discontinuities in the vr*^, 77 and rj' condensates at 

6' = TT (Figures 4,5). 
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Vmin 




Figure 4: The value of the potential at the minimum as a function of 6 for the case of 
degenerate u and d quarks, m = 7MeV, rus = 150 MeV, r/S = 200 MeV. 




Figure 5: The condensates a oc 7r°, /3 oc r^, oc r/' as functions of 9. 



If one moves away from 6 = n one of the two solutions of the minimization equation ( |18|) 
becomes the global minimum, whereas the other solution is a local minimum. By crossing 
9 = IT the global minimum switches between the two branches. This description is valid in 
the vicinity oiO = tx. For | O — ti |> rriu/ms the second (local) minimum disappears and there 
is only one (global) minimum. The situation is similar for a small mass splitting between 
rriu and m^, although Witten's condition (|19D is modified for 6 ^ n. Finally, for large mass 
splittings there is only one minimum of the potential for all values of 6 (see also P]). 
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5 Axion Cosmology 



The results described above have some interesting imphcations for axion cosmology, partic- 
ularly for non-trivial configurations like axion domain walls (ADWs). The internal structure 
of the ADW was investigated previously by Huang and Sikivie 0, who were the first to 
notice its internal 7r° structure. In this section we emphasize the additional 77, 77' structure 
of the wall and the nature of the coupling of an ADW to thermal degrees of freedom below 
the QCD phase transition. 

In axion models the role of the 6 parameter is played by a dynamical axion field a, 
whose potential energy is minimized at 9 very close to zero, thus providing a solution to 
the strong CP problem. The models must exhibit an anomalous Peccei-Quinn symmetry, 
f/(l)pQ, at high energies, which is spontaneously broken at some scale fpq. 

Because of the spontaneous breaking of a U{1) symmetry, the cosmology of these models 
is somewhat involved. At temperatures below /pg, global cosmological strings are formed, 
which persist until the QCD chiral phase transition (unless inflation occurs with reheat 
temperature below /pg). At the QCD phase transition the degeneracy corresponding to 
different values of (a) / fpq ~ 6' is lifted, with the energy minimized at ^ = 2?™ , n E Z in 
the simplest (N=l) axion models. Since the axion field undergoes a 'winding' about the axis 
of the string, after the chiral phase transition the axion strings become connected by domain 
walls in which the axion field interpolates from 9 = Q to 9 = 2tt over a lengthscale m'^^ . 
Our earlier results imply that there are also domain walls in the vr*^, 77 and rj' condensates 
across the ADW. We note that the physical axion and Goldstone fields are obtained by 
diagonalizing the low energy mass matrix, see for example 

The subsequent evolution of a network of axion strings connected by domain walls was 



studied in [|T0|, with the conclusion that the ADWs could persist for many horizon times 
and eventually even affect structure formation. This conclusion relies on the assumption 
that the ADW is extremely weakly coupled to the thermal background particlesp]. In 
this assumption was justified by detailed calculations. However, only the tt*^ component of 
the ADW was considered in 0, and in particular only the derivative couplings of the vr*^ 
to nucleons. The derivative couplings lead to effects which are suppressed by the spatial 
derivative of the wall, which is ~ ma and hence very small. 

There are additional interactions of the ADW with the thermal background which lead 
us to believe that, in fact, the ADW is strongly coupled and can rapidly dissipate excess 
energy. The relevant couplings are due to 

• The internal 77, rj' structure of the wall, which couples strongly to all hadronic thermal 
modes. In particular, the rj' interactions can be estimated from the large- A^^^ effective 



"it also assumes that an ADW does not rapidly destroy itself by self- intersection. 
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Lagrangian, yielding interactions such as n-point 77' vertices which are suppressed at 
most by powers of l/A'^c, but are otherwise of typical strong interaction size. 



• The non-derivative couplings of the 7r° domain wall, which are suppressed by powers 
of the light quark masses (explicit chiral symmetry breaking) but which are still non- 
negligible. For example there are four tt interactions with coefficient ~ Sm//^. 

Due to the interactions mentioned above, the ADW will appear as a potential well of 
depth Aqcd to the thermal hadrons at the QCD phase transition, which we expect to scatter 
with probability of order one. The ADW will contract at a rate 

p/2 /21) 

dt dt ^ ^ ^ 

where / is the size of the ADW, a ~ fpQf-wfTi-K is the surface tension of the domain wall and 
p the hadron energy density. We will focus on the contribution to (^) from pions, because 
they are likely to have the highest energy density after the QCD phase transition. Since the 
quark condensate approaches zero at the chiral phase transition, we expect that the effective 
masses of the pions should be reduced. We therefore estimate the pion energy density to be 
given roughly by the standard relativistic formula for massless {m << T) particles 

P. = "^T' ~ /^^ . (22) 

This should be compared with the energy density of axions which is much smaller: 

Pa ~ 10-V>' • (23) 

The contraction rate of the ADW from interactions with the vr and t]' thermal background 
is therefore very large 

It ''Wo- 

The rate of shrinkage in (|2^) is to be compared with the horizon doubling time ~ Mp/f^. 
Since the Planck scale Mp » fpq, we expect the ADW to contract away in a relatively 
short cosmological timescale, even if the effects of self-intersections are ignored. 
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